Aji and Goldenfeld Reply: In a recent letter wc 
analyzed the critical dynamics of the superconducting 
to normal phase transition in zero magnetic field. We 
explained Monte Carlo (MC) simulation results Q in 
both strong and weak screening limits, where the dy- 
namic critical exponent was found to be zmc ~ 2.7 and 
zmc ~ 1.5 respectively. These results, taken at face 
value were surprising, departing strongly from scaling 
expectations based on model A dynamics that z ~ 2 
H . We showed that the simulations do not measure the 
true dynamic exponent z and that in both the short and 
long range limits, provided the screening length is smaller 
than the system size, the dynamic exponent correctly in- 
ferred from the simulations is z ~ 2, thus removing the 
discrepancy. 

In the preceding Comment Q , the author asserts that 
in our theory, the vorticity has the wrong scaling dimen- 
sion, inconsistent with standard scaling theory and not 
supported by MC data. Second, he asserts, without jus- 
tification, that the identification of MC time with real 
time is correct, and that our proposal that they are not 
the same is an influence of the discreteness of the lattice. 

In fact, the non-equivalence of MC time and real time 
is a genuine effect, but arises due to an implementation 
of a dynamic MC algorithm that does not properly ac- 
count for the scaling of J. Here, we show that our theory 
predicts a scale dependent coupling constant J(L), and 
that the standard finite size scaling results are indeed 
satisfied, contrary to the statement in the Comment. 

In terms of the vorticity, n = V x V0, the free en- 
ergy is, F — (3J2i j ni.njGij[Xo\, where the lattice Green 


function is 


exp[ifc • (fl - fj)] 
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Here J is the coupling constant, Ao is the screening length 
and [3 = 1/fcsT. In the long range case, Ao — > oo, 
the Gij ~ J/l 7 ^ ~ while in the short range case, 


G 


JS(r 


j). The free energy density scales as 


L~ d (F ~ 9, the phase, has a trivial scaling di- 

mension of and n ~ £ -2 . In the long range case, since 
Gij ~ J/£, we get J ~ £, 2 /L 3 where £ is the correlation 
length ((L/£) d ~ L 2d £~ 4 J/£). In the short range case, 
where Gij ~ J/L d , a similar analysis yields J ~ £. In 
the dynamic MC simulations, the scaling of J is not ac- 
counted for and leads to scaling dimensions of n other 
than —2. If one ignored the scaling of J, the above anal- 
ysis would yield n ~ ^~ x and V# ~ £ t 1 ~ x , where x is 
5/2 {(L/£) d ~ L 2d £- 2x /£) in the weak screening limit 
and x = 3/2 in the strong screening limit. The super- 
fluid density is obtained from, F = (3 J drp\V9\ 2 . Given 
the scaling dimension of V# above, and assuming that 
J does not scale, p ~ ((£/£)<* _ pL d^2-2x^ 

which would disagree with standard scaling theory || . 

This discrepancy in the scaling dimensions arises only 
in MC simulations where a time dependent variable is 
measured, such as a two-time correlation function. The 


superfluid density in static MC simulations, for example, 
is obtained from the equal time current-current correla- 
tion function in the phase representation of the Villain 
model. Unlike the vorticity, the phase variable does not 
scale with system size (9 ~ L°). This computation of 
the correlation function reproduces the standard scaling 
form for p. 

The non-equivalence of real time and MC time arises 
from using the same value of J for all lattice sizes L ||, 
and as explained in our letter, is a consequence of the 
fact that in a single MC time step, regardless of L, an 
equal change must be made to the voltage pulse produced 
by the vortex loops and not the vortex loops themselves. 
In dynamic MC simulations, if we scaled J appropriately, 
we would obtain z ~ 2, for strongly screened interactions. 
However, one usually does not know ahead of time what 
the scale dependence of J(L) is. 

In our theory, the scaling form for p is calculated from 
p ~ J(L)L 2x ~ d ~ 2 . In the case of short ranged interac- 
tion and the physically relevant case of weak screening 
(A < L), x = 3/2 and J - L (f - L at T c ), we get 
p ~ L 2 ~ d as expected. In the long range case, where 
Ao is set to infinity, the same argument yields p ~ L 2 ~ d 
since x = 5/2 and J ~ L^ 1 in this case. For the mag- 
netic permeability, including the scaling of J implies, 
p ~ J _1 L d_2;r . With this correction, we obtain the cor- 
rect scaling form, p ~ L d ~ 4 , in both the strong and weak 
screening limit. 

In summary, our analysis is fully consistent with stan- 
dard results, and explains the surprising results of the 
MC simulations. 
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